Study of the pp — ► npiT + reaction at 1.25 GeV with HADES 
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Abstract 

In pp collisions at 1.25 GeV kinetic energy, the HADES collaboration aimed 
at investigating the di-electron production related to A (1232) Dalitz decay 
(A + — ► pe + e~). In order to constrain the models predicting the cross section 
and the production mechanisms of A resonance, the hadronic channels have 
been measured and studied in parallel to the leptonic channels. The analyses 
of pp — ► np7r + and pp — > ppir channels and the comparison to simulations 
are presented in this contribution, in particular the angular distributions being 
sensitive to A production and decay. The accurate acceptance corrections have 
been performed as well, which could be tested in all the phase space region 
thanks to the high statistic data. These analyses result in an overall agreement 
with the one-7r exchange model and previous data. 

1 Introduction 

The High-Acceptance Dielectron Spectrometer (HADES) is a magnetic spectrometer in operation at 
the heavy ion synchrotron facility SIS at GSI Darmstadt. It is designed for di-electron spectroscopy in 
nucleus-nucleus collisions around 1-2 AGeV to study hot and dense nuclear matter. Di-electron pro- 
duction in C+C collisions at 1 and 2 AGeV have been reported on [1, 2]. In the invariant mass region 
0.15 GeV/c 2 < M e+e - < 0.5 GeV/c 2 the measured pair yield shows a strong excess above the contribu- 
tion expected from hadron decays after freeze-out. In this mass range, the r\ Dalitz decay is well under 
control, but two other di-electron sources which play an important role are poorly known: the A (1232) 
resonance Dalitz decay (A + — > pe + e~) and NN bremsstrahlung (AW — » NNe + e~). 
For a better understanding of the contribution from such processes, HADES studied pp and dp inter- 
actions at 1.25 GeV kinetic energy [3,4], just below the r\ production threshold. Due to the small ex- 
pected contribution for the AW bremsstrahlung, the pp reaction at 1.25 GeV is expected to be mostly 
sensitive to the A Dalitz decay [5]. This process is studied by leptonic inclusive (pp —>■ e + e~X) and 
exclusive {pp — ► ppe + e~) channel analyses. At the same time, we measured the hadronic channels 
(pp — ► npir + , pp — > pp-K ) which provide a consistency check with the analysis of leptonic channels as 
well as a precise and independent control of inputs (cross section, angular distribution, etc.) for A Dalitz 
decay study. The pion production is dominated by A excitation at T^i n = 1.25 GeV, therefore, the two 
hadronic channels are correlated by isospin symmetry: a(pp — ► nA ++ ) ~ 3 (pp — > pA + ), leading to 
a(pp — > npi: + ) ~ 5 (pp — > ppn°). 

2 Experiment 

The HADES detector, as shown in Fig.Q} consists of 6 identical sectors covering the full azimuthal range 
and polar angles from 18° to 84° with respect to the beam direction. Each sector contains: A Ring 
Imaging CHerenkov (RICH) detector used for electron identification; two sets of Mini-Drift Chambers 
(MDC) with 4 modules per sector placed in front and behind the magnetic field to determine momenta of 
charged particles; the Time Of Flight detectors (TOF/TOFINO) and the Pre-Shower detector improving 
the electron identification. For reaction time measurement, a START detector is located in front of the 
target. An (e + ,e~) invariant mass resolution at the lo peak of ~ 2.7% and a momentum resolution 
for protons of ~ 3% can be achieved. The first level trigger is obtained by a fast multiplicity signal 
coming from the TOF/TOFINO wall, combined with a reaction signal from the START detector, while 
the second level trigger is made by using the informations from the RICH and Pre-Shower to enrich the 
lepton candidates signals. 

In the April 2006 experiment, we used protons at 1 .25 GeV kinetic energy, with an intensity of about 10 7 
particles/second, and a target filled with liquid hydrogen. The START detector was not used in this run 
because of the too high secondaries produced by the START detector itself. Thus a specific algorithm [6] 
has been developed to calculate the time of flight and then to identify the charged particles. 
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Fig. 1: Schematic layout of the HADES detector. 



3 Simulation 

The two hadronic channels, pp — > npir + and pp — > pp7r°, were simulated at a beam energy of 1.25 
GeV using the PLUTO event generator [7]. As summarized in table [T] the channels involving A and 
N* resonances are simulated with the cross sections taken from [5]. For the A resonance production, 
the mass distribution and angular distribution are taken from Dmitriev's calculation based on one-7r 
exchange model which describes very well the previous experimental data [8]. For the A hadronic 
decay, we implemented the angular differential cross section du/dQ ~ (1 + 1.35 cos 2 #) which is in 
agreement with the anisotropics measured by previous experiments [9, 10] (0 is the 7r polar angle in A 
frame with respect to the direction of momentum transfer from AT to A in the total center of mass). 
Theoretical predictions [1 1] have been used for the angular distribution of N* production since no data 
exist. Simulated events were filtered through HADES acceptance to compare with the data. 



Table 1: Employed hadronic channels for the pp reaction at 1.25 GeV. 



Outgoing channel 


Cross section 


Production processes 


Corresponding cross section 


pp — ► npir + 


19.4 mb 


pp — > nA ++ 


17.0 mb 






pp — > pA + 


1.9 mb 






pp — > pN* 


0.5 mb 


pp — > ppir° 


4.0 mb 


pp — > pA + 


3.8 mb 






pp — > pN* 


0.2 mb 



4 Results 

4.1 pp — >■ np7T + channel 

The pp — > npir + channel is studied using a reconstruction of the undetected neutron. The reaction 
was selected first by the charged particle identification based on momentum and reconstructed time of 
flight [12], then a (p, n + ) missing mass cut was imposed around the neutron mass. This cut efficiently 
suppresses the background coming from misidentified protons and two-7r contributions. In the Dalitz plot 
(Fig f2l), one can clearly see the A ++ signal located around M? nv (p, n + ) = 1.5 (GeV/c 2 ) 2 corresponding 
to the squared mass of A ++ (1232), while the A + signal located around M^ nv (n, ir + ) = 1.5 (GeV/c 2 ) 2 is 
less pronounced. The spot appearing in the Dalitz plot for (n, 7r + ) and (p, tt + ) invariant masses squared 
around 2 (GeV/c 2 ) 2 is due to the pn Final State Interaction (FSI). This is confirmed by a simulation 
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Fig. 2: Dalitz plot of the pp 
(preliminary HADES data). 
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where pn FSI was implemented using the Jost function formalism [13]. 

Fig. [3] exhibits the projection on the (p, tt + ) (Left) and (n,7r + ) (Right) invariant masses and a 
comparison to simulation. The data and simulation are both normalized to the total pp elastic cross 
section which is also measured in this experiment. Error bars include statistical and systematic errors 
due to event selection (5%) and correction of trigger condition which has been applied for some of the 
bins. The uncertainty on the normalization to pp elastic scattering is also considered as a source of 
systematic error (6%) but not included in the error bars here. The prominent peak of M- mv (p, tt + ) around 
1.23 GeV/c 2 confirms that most of the n + are produced via A ++ decay, which is consistent with the 
resonance model. N* contribution seems also reasonable and the invariant mass distributions are rather 
well reproduced by the pp — > nA ++ and pp — > pA + simulations. 
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Fig. 3: Left: (p, ir + ) and Right: (n, ir + ) invariant mass distributions compared to Pluto simulation: total contri- 
bution (red), A ++ (blue), A + (pink) and N* (green). Both the data and simulations are normalized to the total pp 
elastic cross section. 



We also looked at the neutron angular distribution in the center of mass system (Fig. HJ) which 
is sensitive to the angular distribution of A resonance production since pp —>■ nA ++ is the dominant 
process. The comparison shows a good agreement between data (black points) and simulation (red solid 
line) for forward and backward neutron angles, however, an excess around cos# n = is found (our data 
are a factor 2 above simulation). By adding a non-resonant ir production contribution with a = 0.5 mb 
the simulation (red dashed line) got closer to the data points around 9 n = 0, even though, it can not fully 
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explain the discrepancy. On the other hand, the cross section assumed here for the non-resonant contri- 
bution seems reasonable and can not be larger since it is also constrained by invariant mass distributions. 
But one should mention that the majority of statistics is located at forward and backward neutron angles 
where the yield is in agreement with the simulation and that the overall discrepancy is less than 5% of 
the total cross section. To minimize the model dependence, a two dimensional acceptance correction has 
been performed using the resonance model. For each (cos# n , M- mv (p, 7r + )) bin, the correction factors are 
calculated as the ratio of events from simulation in full phase space and in geometrical HADES accep- 
tance, then they are applied to correct the data. The width of the bins is adjusted to optimize the precision 
of the correction. The systematical error introduced by the acceptance correction procedure is estimated 
to 5%. We finally obtain the cross section a(pp — ► npir + ) = 20.4 ± 1.9 mb, in good agreement with 
previous data [5]. 
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Fig. 4: Left: Angular distribution of neutron in the pp center of mass system. Data (black points) compared to 
simulations. (Pluto simulation (red solid line) including A ++ (blue), A + (pink) and N* (green); improved Pluto 
simulation (red dashed line) with additional non-resonant contribution (orange dashed line)). Right: Angular 
distribution of neutron in center of mass system after acceptance correction. Data (black points) compared to Pluto 
simulation (red line). Both data and simulations were normalized to the total pp elastic cross section. 



0.8 




Fig. 5: Left: Angular distribution of ir + in (p, ir + ) center of mass system. Experimental data are presented as black 
points and compared to simulations with different A decay anisotropy parameters (see text). Pluto simulation: 
B = 1.35 (red solid line), pure one-7r exchange model: B = 3 (red dot-dashed line) and isotropic: B = (red 
dashed line). Right: Angular distribution of tt + in (p, ir + ) reference frame after acceptance correction. Data (black 
points) are compared to Pluto simulation (red line). 



The decay angular distribution of the A resonance has been studied in this channel as well by 
looking at the angular distribution of 7r + in the (p, n + ) center of mass system. In the simulation, the 
angular distribution of 7r in the A resonance frame according to da/dQ ~ (1 +B cos 2 8) is assumed. In 
order to see the sensitivity of the 7r angular distribution to the anisotropy of the A decay, and in addition to 
the decay distribution mentioned above in Pluto simulation (where B = 1.35), the distributions expected 
for pure one-7r exchange model (where B = 3) and isotropic decay (where B = 0) were simulated 
and filtered through the HADES acceptance. The comparison to the data (left panel of Fig. [5]> shows 
the great sensitivity to the anisotropy of the A and the good agreement with anisotropics measured in 
previous experiments [9, 10]. The acceptance correction is done using the same procedure as for the 
neutron angular distribution mentioned above. The right panel of Fig. [5] shows the angular distribution 
of 7r + in the (p, ir + ) reference frame extrapolated to full solid angle, which fits nicely with the simulation. 

4.2 pp — > pp-K channel 

The ppir reaction was extracted from (p, p) missing mass spectrum by subttacting a background under 
the 7T° peak. The left panel of Fig. |6]presents the Dalitz plot of this channel. The prominent A+ signal 
(right panel in Fig. |6]) is in good agreement with the event distribution obtained from the resonance model 
with cross section in Table [T] This analysis results in an exclusive ir° cross section extrapolated to 4tt: 
a(pp — > ppit ) = 4.05 ± 0.37 mb, in agreement with previous data. The estimated a(pp — > npTT + ) ~ 
5 (pp — > ppn ) is well reproduced [5]. 




Mf nv (p,7t°) [(MeV/c 2 ) 2 ] M f™P"° [(MeV/c 2 ) 2 ] 

Fig. 6: Left: Dalitz plot of the pp — > ppn° reaction: (p, ir°) invariant mass squared distributions with preliminary 
HADES data. Right: Projection of the Dalitz plot onto the (p, it ) invariant mass squared (black points) compared 
to the simulation (red) assuming the A + (green) and the N* (blue) as intermediate resonances. Both disUibutions 
are normalized to the same number of pp elastic events. 



5 Conclusion and outlook 

In summary, the preliminary HADES results for hadronic channels pp —> npir + and pp — > ppir at 
1.25 GeV have been presented. Both measured it production cross sections agree with the previous 
data and the resonance model within error bars. The invariant masses and angular distributions show an 
overall agreement with the one-7r exchange model although a deviation is observed at large production 
angles. Consistent with previous data, the A resonance decay is clearly anisotropic. With the large 
statistics achieved, the spectra could be extrapolated to full phase space using 2-dimensional acceptance 
correction with a minimized model dependence. This provides a useful test of the resonance model used 
for the analysis of the di-lepton channels and is especially important for the ongoing study of A Dalitz 
decay, using (pe + e~) events. 
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